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Abstract 

The decay B*^ D^D^ offers an interesting probe of CP violation, but requires 
control of penguin effects, which can be done through B'^ D^D~ by means of the 
f7-spin flavour symmetry of strong interactions. Recently, the Belle collaboration 
reported indications of large CP violation in the i?^ decay, which were, however, 
not confirmed by BaBar, and first signals of the B^ channel were observed at 
the Tevatron. In view of these developments and the quickly approaching start 
of the LHC, we explore the allowed region in observable space for CP violation 
in B^ D^D~, perform theoretical estimates of the relevant hadronic penguin 
parameters and observables, and address questions both about the most promising 
strategies for the extraction of CP-violating phases and about the interplay with 
other measurements of CP violation and the search for new physics. As far as the 
latter aspect is concerned, we point out that the B^ D^D~ system provides 
a setting for the determination of the B^-B^ mixing phases {q G {d,s}) that is 
complementary to the conventional -B^ J/ipKs and B^ J/ip'P modes with 
respect to possible new-physics effects in the electroweak penguin sector. 
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Figure 1: Tree and penguin topologies contributing to the [/-spin-related D~^D~ 
and Bg DfD~ decays. 

1 Introduction 



Decays of B mesons are subject of intensive investigations. Thanks to the interplay 
between a lot of theoretical work and the data from the e~^e~ B factories with their 
detectors BaBar (SLAG) and Belle (KEK) as well as from the Tevatron (FNAL), valuable 
insights into CP violation could be obtained during the recent years. In the Standard 
Model (SM), this phenomenon is closely related to the Cabibbo-Kobayashi-Maskawa 
(CKM) matrix [1], and can be characterized by the unitarity triangle (UT) with its 
three angles a, P and 7. This adventure will soon be continued at the LHC (CERN), 
with its dedicated 5-decay experiment LHCb. 

In the 5-physics landscape, an interesting probe of CP violation is also offered by 
B^ —>■ D^D^. As can be seen in Fig. [U this decay originates from b — > ccd quark- 
level processes, and receives contributions both from a colour-allowed tree-diagram-like 
topology and from penguin diagrams. In analogy to the prominent B^ tt+tt" decay, 
the latter contributions lead to complications of the theoretical interpretation of the 
CP-violating observables. However, the penguin effects can fortunately be controlled 
by means of the 5° — > DfD~ channel [2], which is related to B^ —>■ D^D~ through 
an interchange of all down and strange quarks, as can also be seen in Fig. [H Because 
of this feature, the ?7-spin flavour symmetry of strong interactions allows us to derive 
relations between non-perturbative hadronic parametersjl] so that the measurement of 
CP violation in B^ — »• D~^D~ can be converted into CP- violating weak phases. In 
comparison with conventional flavour-symmetry strategies [3], the advantage of the U- 
spin method is that no additional dynamical assumptions are needed, and that also 
electroweak (EW) penguin contributions are automatically included. 

The key observables are the CP-averaged branching ratios as well as the direct and 
mixing- induced CP asymmetries Aqp^B^ — > D^D~) and A^^{Bd D~^D~), respec- 
tively, which enter the following time-dependent rate asymmetry [4]: 



Acp{Bd{t)^D+D- 



TjB^t) ^ D^D-) -r{B%t) ^ D^D- 
r(50(t) ^ D+D-) + r(50(t) ^ D+D- 



^The C/-spin flavour symmetry connects strange and down quarks in the same way through SU{2) 
transformations as the isopsin symmetry connects the up and down quarks. 
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= A%l.{Bd ^ D+D-) cos(AMrft) + A^^Bd ^ D+D-) sin(AM,t), (1) 

where AM^ is the mass difference of the Bd mass eigenstates. The Belle collaboration 
has recently reported evidence for CP violation in D^D~, which could not be 

confirmed by BaBar. The current status reads as follows: 

Jdir/R / +0.11 ±0.22 ±0.07 (BaBar [5]) , . 

•^cp[^d ^ ^ ^ J - I _Q ^ Q 23 ± 0.06 (Belle [6]) ^ ^ 

jdir/R _ n+n-\ _/ +0.54 ±0.34 ±0.06 (BaBar [5]) , . 

^cp[^d^iy u ^ - \+1.13± 0.37 ±0.09 (Belle [6]); ^"^^ 

the Heavy Flavour Averaging Group (HFAG) gives the following averages [7]: 

A%l,{Bd^ D+D-) = -0.^7 ±Q.l7, ^^^"(5d ^ = 0.75 ± 0.26, (4) 

which have to be taken with great care in view of the inconsistency between the BaBar 
and Belle measurements. Concerning the CP-averaged branching ratio, we have 

RR^R ^n+n-N_/(2.8±0.4±0.5)xlO-^ (BaBar [8]) , 
^ ^d^d) I ^197 ± 0.20 ± 0.20) X 10-4 (Belle [6]), ^ ' 

yielding the average of BR(Sd DjD^) = (2.11 ± 0.26) x 10"^ Thanks to the 
updated Belle result, this number is now about 1.6 a lower than the HFAG value of 
BR{Bd DjD-) = (3.0±0.5) x 10"^ [7]. The CDF collaboration has recently observed 
the first signals of the B^ D^D^ decay [9], which correspond to the CP-averaged 
branching ratio 

BR(5, ^ D+D-) = (1.09 ± 0.27 ± 0.47)%. (6) 

Performing a run on the T(55') resonance, also the Belle collaboration has recently 
obtained an upper bound of 6.7% (90% C.L.) for this branching ratio [10]. Moreover, 
the DO collaboration has performed a first analysis of the combined Bg Dl Dl 
branching ratio, with the result of BR(5, ^ D'f'^o'f^) = (3.9lJ:?l}|)% [11]. 

Although the current experimental picture is still in an early stage, it raises several 
questions, which are further motivated by the quickly approaching start of the LHC: 

• What is the allowed SM region for the CP violation in B^ — * D^D^l 

• What are the most promising strategies for the extraction of weak phases? 

• What is the interplay with other measurements of CP violation and the search for 
new physics (NP)? 

These items are the central target of this paper. It is organized as follows: in Section [21 
we explore the parameter space of the CP-violating B^ D^D~ asymmetries, taking 
also the constraints from and similar modes into account, and perform 

a theoretical estimate of the corresponding observables in Section [31 In Section [H we 
discuss the extraction of CP-violating phases from the B^ D^D~ and 5° D^D~ 
decays, while the interplay with other CP probes is discussed in Section [5l Finally, we 
summarize our conclusions in Section [6l 
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Figure 2: The situation in the A^^{Bd D+D-)-A'^^{Bd D+D-) plane. Left 
paneh contours following from the general SM parametrization; right panel: constraints 
following from a measurement of the quantity H. 



2 CP Violation in D+D' 
2.1 Standard Model Expressions 

In the SM, we may write the D^D^ decay amplitude as follows [2]: 

^ D+D-) = -XA [1 - ae'^e'^] , (7) 

where A is the well-known Wolfenstein parameter of the CKM matrix [12], A denotes a 
CP-conserving strong amplitude that is governed by the tree contributions, while the CP- 
consering hadronic parameter ae*^ measures - sloppily speaking - the ratio of penguin 
to tree amplitudes. Applying the well-known formalism to calculate the CP-violating 
observables that are provided by the time-dependent rate asymmetry in ([T]), we obtain 
the following expressions: 

A^^iB, D-D-) = '"""^""^ ^ (8) 

' 1 -2acos^cos7 + a2 ^ ' 

yimix/R n+-n-\ sin 0^ - 2a cos 6 sin(0d + 7) + sin(0rf + 27) 

1 — 2a cos cos 7 + a"' 

where (pd denotes the B^-B^ mixing phase, which takes the value of 2(3 in the SM. 
This quantity has been measured at the B factories with the help of the "golden" decay 
B^ J/ipK^ and similar modes, including Bd J/ipK* and Bd D*D*Ks channels 
to resolve a twofold ambiguity, as follows [7]: 

(Pd = (42.6 ± 2)°. (10) 

Concerning the angle 7, the SM fits of the UT obtained by the UTfit and CKMfitter 
collaborations [13, 14] yield 7 = (64.6 ± 4.2)° and 7 = (59.0l|?)°, respectively. A 
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recent analysis of the f/-spin-related — >■ tt+tt^ and Bg K'^K^ transitions finds 
7 = (66.615 q^3'q)° [15], in excellent agreement with these fits. A similar picture emerges 
also from other recent 7 determinations from B -kti^hK decays [16,17]. Thanks to 
the LHCb experiment [18], our knowledge of 7 will soon improve dramatically, also since 
very accurate "reference" determinations through pure tree decays will become available. 
In the limit of a — > 0, expression ([9]) would allow a straighforward extraction of siuf^^. 
However, these penguin effects cannot simply be neglected and require further work. 

For the following discussion, we shall assume 7 = 65° and (pd = 42.6°. Using ([8]) and 
(E]), we can then calculate contours in the A^^{Bd D+D-yA^^iBd D+D-) plane 
for given values of a and 6', which are theoretically exact in the SM. The resulting picture 
is shown in the left panel of Fig. [2] (for its B^ n'^ir' counterpart, see Ref. [19]). There 
we have also included the experimental BaBar and Belle results, as well as the HFAG 
average; the dot-dashed circle defines the outer boundary in this observable space that 
follows from the general relation 

[A^UBd ^ D+D-)Y+ [A'Sl^iBd ^ D+D-)f < 1. (11) 

Fig. [2] shows that the Belle result lies outside of the physical region, in contrast to the 
BaBar measurement and the HFAG average. The contours of that figure allow us to 
read off the corresponding values of a and 6 straightforwardly. 



2.2 Constraints from D+D' 

s s s 

We now go one step further by using the information that is offered by the B'^ 
decay. In analogy to ([7]), its SM amplitude can be written as follows: 



DtD: 





1 + ea'e'^'e'^ 







where 



1- A2 



0.05. 



Following Ref. [2], we introduce 



H = - 

e 



A' 



A 



where 



Mb, <!>{MdJMb„MdJMb,) tb, 

1 — 2a cos 6 cos 7 + 
1 + 2ea' cos 9' cos 7 + e^a'^ ' 

Hx,y) = ^[l-(x + yy] [1-ix 



BR{Bd 



BR(5, 



D+D 



D+D: 



DtD: 



(12) 



(13) 



(14) 
(15) 



is the well-known B — >■ PP phase-space function, and the tb, , are the B^^s lifetimes. 
Applying the ?7-spin flavour symmetry, we obtain the relations 



6' = 9. 



(16) 
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Thanks to the e suppression in (fT4|) . the impact of f/-spin-breaking corrections to (fT6l) 
is marginal for H. In the case of |^Y^|, ratios of fZ-spin-breaking decay constants and 
form factors enter. If we apply the "factorization" approximation, we obtain 



A! 



A 



_ {Mb, - MpJ ^MbMj,, {ws + 1) fj,^ Uws) 
,,,, {Mb, - Md,) ^Mb.Md, {wa + l) fo, U^a) 



(17) 



where the restrictions form the heavy-quark effective theory for the Bg Dg form 
factors have been taken into account by introducing appropriate Isgur-Wise functions 
^q{wq) with Wq = Mb^/ {2Md^) [20]. Studies of the light-quark dependence of the Isgur- 
Wise function were performed within heavy-meson chiral perturbation theory, indicating 
an enhancement of ^s/^d at the level of 5% [21]. Applying the same formalism to foJfD, 
gives values at the 1.2 level [22], which is in accordance with the recent measurement by 
the CLEO collaboration [23]: 

^ = 1.23 ±0.11 ±0.04, (18) 

as well as with lattice QCD calculations, as summarized in Ref. [24]. Using heavy-meson 
chiral perturbation theory and the 1/A^c expansion, non-factorizable S'[/(3)-breaking 
corrections were found at the level of a few percent in Ref. [25]. The CDF result in ([6]) 
and the average of ([5]) yield then, with the CLEO measurement in f|T8l) . the following 
numbers: 

i/ = 0.59 ±0.31 (0.84 ±0.45), (19) 

where we have added the errors in quadrature, and have also given the result correspond- 
ing to the HFAG value of BR(i?d — > D^D~) in parentheses. The general expression for 
H in f|T^ implies a lower bound [26], which is given by 

H> [1 - 2e cos^ 7 ± 0{e^)\ sin^ 7 0.81. (20) 



Consequently, the rather low central value of (]T9l) . which is essentially due to the new 
Belle result [6], is disfavoured by the experimental information on 7. 

If we replace the s spectator quark of the 5^ DfDj decay through a d quark, 
we obtain the DfD^ process. Whereas the B^^s) ~^ ^d{s)^d{s) system receives 

contributions from tree and penguin as well as exchange {E) and penguin annihilation 
{PA) topologies (the latter are not shown in Fig. [T]), the 5^ — > DfD~ channel and its 
[/-spin partner B^ — ^ D^D^ receive only tree and penguin contributions. Consequently, 
if we use the SU{2)) flavour symmetry and assume that the exchange and penguin anni- 
hilation topologies play a minor role, we may replace 5° — > D^D~ in the determination 
of H through 5° D+D' [27] I Expression ([HD is then modified as follows: 
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^MnjMB,,MnjMBj 
mMBjMB,,MnjMB,)\ 



BR{Bd D+D- 
\BR{Ba^DfD^)\ 



(21) 



^This is analogous to the replacement of i?" K^K through ^ tt [28]. 
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The importance of the E + PA amphtude can actually be probed through the [/-spin re- 
lated -Brf(s) D'^f^^^D~^^^ decays. The current experimental situation can be summarized 
as follows: 

RRfR _^n±n^^_ 1(6.4 ±1.3 ±1.0) X 10-3 (BaBar[29]) 

ati[iJd -^^s^ J - I 5 ± 0.2 ± 0.8 ± 0.8) x lO^^ (Belle [30]), ^ 

yielding the average of BR{Bd DfD^) = (7.1 ± 0.9) x 10"^; Belle reported also the 
upper limit of BR{Bd D+D;) < 3.6 x 10"^ (90% C.L.) [30]. Expression (EID gives 
then 

i/ = 0.85 ±0.19 (1.22 ±0.31), (23) 

where the notation is as in ( |T9l) . 

Let us now investigate the constraints on the A^^{Bd D+D-yA^^{Bd D+D') 
plane that follow from H . If we use ( !T4|) with (fT6|) . we may eliminate the strong phase 
^ in ([8]) and ^ with the help of 

cosg= ^-H + iX-e'H)a^ sin^ = ±Vl - cos^ (24) 
2a(l±ei/)cos7 ' ^ ^ ) 

If we then keep a as a free parameter, we arrive at the situation shown in the right panel 
of Fig. [2], where the dashed line separates the regions with cos^ > and cos^ < 0. 
In the factorization approximation, we expect a negative value of cos 9. Although non- 
factorizable effects could generate a large value of 6*, we do not expect that cos^ changes 
its sign. This feature is in fact observed for other non-leptonic 5-meson decays, such as 
the B^i Tr+TT^, B^ 7t~K^ system [15]. With 7 ~ 65°, which corresponds to cos 7 > 
0, the expression in (IT^ implies then H > 1. In the right panel of Fig. [2], this leaves 
us with the banana-shaped region in the AQp{Bd D'^ D~)-AQp{Bd D^D~) plane. 
Interestingly, the central value of the HFAG average falls well into this region, whereas 
the central value of the BaBar result would require a positive value of cos9. Although 
the current errors are too large to draw definite conclusions, this exercise illustrates the 
usefulness of the plots in observable space to monitor the experimental picture. Since 
the Bg input for the determination of H is just the CP-averaged Bg — * D^Dj branching 
ratio, this measurement would also be interesting for an e+e~ (super-)i? factory operating 
at the T(5S') resonance [10,31]. 



3 Theoretical Estimates 

In order to analyze the B^ — > D^D~ decay theoretically, we "integrate out" the heavy 
degrees of freedom, i.e. the W boson and top quark in Fig. [H and use an appropriate 
low-energy effective Hamiltonian, which takes the following form [32]: 



0/ Gy 



10 



k=l 



k=3 



(25) 



Here the Xf^ = VjaV*^ denote CKM factors, and Qf {3 G {% c}) are the usual 
current-current operators, Qf, . . . ,Qq and Q^, . . . , Q^q denote the QCD and EW penguin 
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Figure 3: Theoretical estimates of the hadronic parameter ae*^ (left panel), and the 
Bd{s) ^d{s)^d{s) observables (right panel) for 7 = 65°, (pa = 42.6° and Rh = 0.45. 



operators, respectively, and /i = 0{mb) is a renormalization scale. If we apply the 
Bander-Silverman-Soni mechanism [33] as well as the formalism developed in Ref. [34], 
we obtain the following estimate: 



ae^^ ^ Rh 



At: + At + Ac_ 

where Rb oc iKib/Kfel is the corresponding side of the UT, and 

At = lc^ + C2 
1 



(26) 



(27) 



At = - {Cs + Cg + XD {C, + Cr)] + + Cio + XD {Ce + Cs) (28) 



Aj 



as 

97T 
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G{mj, k, nib) 



C2 + l-{3C, + C2) 

3 a., 



;i+xd) 



(29) 



with j G {u,c}. The Ck refer to fi = rrib and denote the next-to-leading order scheme- 
independent Wilson coefficient functions introduced in Ref. [35]. The quantity 



Xd 



2Ml 



{rric + md){mb - rric) 



(30) 



is due to the use of the equations of motion for the quark fields, whereas the function 
G{mj, k, TUb) originates from the one-loop penguin matrix elements of the current-current 
operators Q''i2 with internal j quarks. It is given by 



G{mj, k,mb) = — 4: / dxx(l — x)ln 



— k"^ x{l — x) 



mt 



(31) 



where mj is the j-quark mass and k denotes some average four-momentum of the virtual 
gluons and photons appearing in the penguin diagrams [34]. In Fig. [3l we show the 
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Figure 4: Illustration of the contours in the 7-a plane for the central values of the 
CP-violating Bd D^D^ asymmetries in (jlj) and various values of the ratio H of the 
CP-averaged Bd D^D^ , Bs DfDj branching ratios. 

corresponding results, keeping fc^ as a free parameters. The sensitivity on k"^ is moderate, 
and in the case of H and the mixing-induced CP asymmetry even small. It should be 
emphasized that these results, with a ~ 0.08 and 6 ~ 205° yielding the observables 

H ~ 1.07, A'^^iBd D+D-) 5% and A'^^{Bd D+D-) ~ 76%, can only be 

considered as estimates. A similar analysis was also performed in Ref. [36]; however, in 
Eq. (12) of that paper, a factor of ^ is missing in front of C3, and 10/3 should read 10/9. 

It is instructive to compare (1261) with the corresponding expression for the penguin- 
to-tree ratio d of the B^ vr+vr" decay in Ref. [28]. We observe that a is suppressed 
with respect to d by a factor of ~ 0.2. The value of ~ 0.4, as determined from 
the [/-spin analysis of the Bd tt+tt^, Bg K^K~ system [15], points therefore also 
towards a ~ 0.08. However, the detailed dynamics of these decays is of course very 
different, so that values of a at the 20% level cannot be excluded. 

4 Extractions of CP- Violating Phases 
4.1 Extraction of 7 

As was pointed out in Ref. [2] , we may combine ([8]) with ([9]) to eliminate the strong phase 
^, which allows us to calculate a as a function of 7. To this end, the B^-B^ mixing phase 
(pd is needed as an input. The corresponding contour relies only on the SM structure of 
the B^ D^D~ decay amplitude and is theoretically clean. A second curve of this kind 
can be fixed through AQp{Bd D^D^) and H with the help of the f/-spin relations in 
(ITB]) . The advantage of the combination of these observables is that they both depend on 
cos^^'- \ Because of the e suppression of the a' terms in f|T^ . [/-spin-breaking corrections 
to this relation have actually a very small impact, so that the major non-factorizable U- 
spin-breaking effects enter through the determination of H. In Fig. HI we illustrate this 
strategy for the central values of the averages in (jlj) and different values of H . We see 
that H = 1.15 would give a value of 7 = 63° with a = 0.25 (and 6 = 249°). On the other 
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Figure 5: Determination of the hadronic parameter ae*^ for given values of H and 

hand, H = 1.05 yields 7 = 89° with a = 0.22 and 9 = 244°, whereas H = 1.25 results 
in 7 = 42°, a = 0.35 and 6 = 257°. Consequently, since a variation oi H = 1.15 ± 0.10 
gives the large range of 7 = (68^21)°, the situation would not be favourable for the 
determination of this UT angle. However, the hadronic parameter a = 0.251o'q3 and 
- in particular the strong phase 9 = (249l5)° ^ could be well determined, but are of 
less interest. In the case of the f/-spin- related 5^ tt+vt^, Bg K^K~ decays, the 
current data result in a complementary situation, with a very favourable situation for 
the extraction of 7, and a less fortunate picture for the corresponding strong phase [15]. 
It will be interesting to follow the future evolution of the Bd{s) ^d{s)^d{s) data. 

4.2 Extraction of the B^B^ Mixing Phase 

An alternative avenue for extracting information from the CP-violating asymmetries of 
the B^ — s> D^D^ decay arises if we use 7 as an input. By the time accurate measurements 
of these CP asymmetries will become available we will also have a clear picture of this 
UT angle thanks to the precision measurements that can be performed at LHCb [18]. 
For the following analysis, we assume a value of 7 = 65° (see the remarks after (fTO!) ). 
If the penguin effects could be neglected, the following simple situation would arise: 

(sin 2f3)D+D- = ^cp (^rf ^ D+D-) sin 0^ ^= sin 2(3. (32) 

The goal of the following discussion is to include the penguin effects in the determination 
of sin (pd- To this end, we first determine a through the combination of Aq^IB^ —>■ D~^D~) 
and H by means of the f/-spin relation ( |T6l) . which yields 

a = \/b-Vb'^-l, (33) 

where 

6iV = 2 [(1 + ei7) sin 7 cos 7] V (i/ - 1) (1 - e'^H) sin^ 7 

-e [(1 + e)HA'^l>iBd ^ D+D-) cos 7]' (34) 
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Figure 6: Correlation between A^^{Bd — > D^D^) and sin^^ for given values of H and 
various values of A%l.{Bd D+D-): (solid), ±0.1 (dotted), ±0.2 (dashed), ±0.3 



(dot-dashed). The shaded region corresponds to the experimental value of (sin 2/3) 



cN= [(i7- l)sin7]V [{1 + e)HA'^^{Bd D+D-) cos 



(35) 



with 



AT = [(1 - e^H) sin7] ± [e(l ± e)HA'^'p{Bd D+D-) C0S7] . (36) 

In f l33|) . the sign in front of the inner square root could, in principle, be positive or 
negative. However, since the large values of a corresponding to the ± sign are completely 
unrealistic, we have already written the — sign. The strong phase 6 follows then from 

1 - ^ + (1 - e'H)a' ^g^^ 



sin^ 



cos 9 

2(1 ± eiJ)a cos 7 
(l + e)(l±ea2 



HAcp 



(Bd ^ D+D 



(38) 



2(1 ± eiJ)asin7_ 

In these expressions, the impact of the e terms is tiny, but we have kept them for 
completeness. In Fig. |5l we show the resulting picture of the hadronic parameter ae*^ 
in the complex plane for various values of H and AQp{Bd D^D^), which should be 
compared with theoretical estimate shown in the left panel of Fig. [31 

If we now use again fll6p and eliminate cos 6' in through (I37p . we obtain 



A sin (f)d + B cos ipd = C, 



where 



A = [H - 2a^ sin^ + eH {1 + (l - 2 sin^ 7 ± e) a^}] C0S7 
B = [H - 1 + cos2-f + eH {1 + cos 27 ± e) a^] sin 7 
C = (1 + e)(l ± ea^)HA'S^{Bd D+ D') cos 7, 
with a given in (l33ll . Finally, sin^^ can be determined as follows: 

AC - BVA^ + B^-C^ 



sm (pd 



A'^ + B^ 
10 



(39) 

(40) 
(41) 
(42) 

(43) 



Here we have chosen the sign in front of the square root such that we obtain a positive 
value of cos0rf, in agreement with the 5-factory data for the CP-violating effects in 
the Bd —>■ J/tpK* and Bd —>■ D*D*Ks channels [7]. In Fig. [6l we show the resulting 
correlation between the mixing-induced CP violation in B'^ —>■ D^D^ and sin^^ for 
various values of H and A^p^Bd D^D~), which correspond to the situation shown 
in Fig. [5l These curves allow us straightforwardly to include the penguin effects in 
the determination of the B^-B^ mixing phase form the CP-violating B^ —>■ D^D~ 
observables. 



4.3 Extraction of the B^-B^ Mixing Phase 

Let us now turn to the CP-violating rate asymmetry of the B^ 
is defined in analogy to ([1]), and takes the form 



D^D^ decay, which 



A^ABs^DtD- 



cos{AM,t) + A^'^{B, 



DtDl 



sin(AMst) 



cosh(Ar,t/2) - AaABs D+D-) sinh(Ar,t/2) 



(44) 



where AF^ = F^ ^ — F^ is the difference of the decay widths F^'' and F^ of the "heavy" 
and "light" mass eigenstates of the Bg system, respectively. The mass difference AM^ 
was recently measured at the Tevatron [37,38], with a value that is consistent with the 
SM expectation. On the other hand, this result still allows for large CP-violating NP 
contributions to B^-B^ mixing (see, for instance, Refs. [39,40]). In this case, the mixing 
phase (ps would take a sizeable value, and would manifest itself also through significant 
mixing- induced CP violation in B^ D^D' at LHCb. In the SM, we have on the other 
hand a tiny phase of 0^ = —2X^7] —2°, where i] is another Wolfenstein parameter. 
Using the formalism discussed in Ref. [4], (fT2l) yields 



A!^p{Bs 



v4.Ar(-Bs 



D+D: 



d+d: 



2ea' sin 6' sin 7 



1 + 2ea' cos 6' cos 7 + e^a'^ 
sin (ps + 2ea' cos 9' sin(0s + 7) + e^a'^ sin( 



+ 27) 



1 -|- 2ea' cos 6' cos 7 + e^a'^ 
cos (j)s + 2ea' cos 6' cos(0s + 7) + e^a'^ cos( 



+ 27) 



1 + 2ea' cos 0' cos 7 -|- e^a'^ 
and f|T6|) implies the following ?7-spin relation [2]: 

,(5. ^ DtD-) 



(45) 
(46) 
(47) 



A'^^{Bd 



D+D- 



-eH. 



(48) 



Thanks to the suppression through the e parameter in fH6l) . the penguin effects are 
significantly smaller than in the case of B^ D^D~. Nevertheless, since we are aiming 
at precision measurements, it is important to be able to control them. Since we may 
determine the penguin parameters a and 9 as we have discussed above, the [/-spin 
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relations in (fT6|) allow us to include the penguin effects also in the determination of 0^- 
It is instructive to perform an expansion in powers of ea', which yields 

sin0, = A^^{Bs DtDj) T 2ea'cos^'sin7-^l - A^'^{B, D+D- f + C((ea')'), 

(49) 

where ^ refers to sgn(cos0s) = ±1. For strategies to determine this sign, which is 
positive in the SM, see Refs. [15,41]. Using fl37l) . the relevant hadronic parameter can 
straightforwardly be fixed: 

2ea'cos^'sin7 = (1 - i/) tan7 + ^(a^). (50) 

Let us finally have a closer look at the observable 

^Ar(5, DjD;) = -COS0, + 2ea'cos0'sin7sin0, + C((ea')^), (51) 

which can be extracted from the following "untagged" rate: 

{T{B,{t) ^ DtD:)) = r(i?°(t) ^ DjD;) + Vmt) DtD-) 

oc e-r^* [e+^^^*/2^L(i?. ^ DjD;) + e'^^^^'/^RniB, ^ D+D-)] . (52) 

Here Tg denotes the average of the decay widths of the "heavy" and "light" mass eigen- 
states of the Bg system, and 

i?L(5, D+D-) = 1 - AAr{Bs D+D-) = 1 + cos 0, + C(ea') 2, (53) 

Rn{Bg ^ D+D-) = 1 + AAriBg ^ D+D;) = 1 - cos 0, + O(ea') ~ 0. (54) 

As far as a practical measurement of fl52|) is concerned, most of the data come from short 
times with ATst ^ 1. We may hence expand in this parameter, which yields 

{T{B,{t)^DjD;))^e-^^' 

Moreover, if the two-exponential form of (|52l) is fitted to a single exponential, the corre- 
sponding decay width satisfies the following relation [41]: 

AT 

^DtDT = r. + AAviBs ^ DjD;)^ + o((Ar,)Vr,). (56) 

Using flavour-specific B^ decays, a similar analysis allows the extraction of up to 
corrections of C((AF,/F,)2) [41]. In the presence of NP, AF^ is modified as follows [42]: 

AF, = AFS^^cos0„ (57) 

where AF^^/F^ is negative for the definition given above, and calculated at the 15% 
level [43]. Consequently, (156!) actually probes 

|AF^^| 

F^+ jj- -Ts = [cos^ (j)s - ea cos 9' sin(20s)] — ^ h . . . , (58) 

thereby complementing other determinations of the width difference of the Bs system, 
such as from the ?7-spin- related Bg — >■ K~^K~, Bd tt+tt" decays [15]. 



Aat{.Bs 



DtD: 



+0{{AVstf) . (55) 
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5 Interplay with Other Probes of CP Violation 



As we have seen in the previous section, the ?7-spin-related Bq ^q^q decays offer 
an interesting tool for the extraction of the B^-Bg mixing phases (g G {d,s}). Since 
the "golden" decay B^ — > J/ipK^ and similar channels allow already a very impressive 
determination of 0^, as can be seen in (11 01) . this may not look as too exciting. However, 
this is actually not the case. In fact, the current value of (fTOj) is on the lower side, and 
the interplay with the UT side Rh oc |\4fe/Kf,| leads to some tension in the CKM fits 
[7, 13, 14], which receives increasing attention in the 5-physics community. If this effect 
is attributed to NP, the standard interpretation is through CP-violating contributions 
to 5^-52 mixing, with a NP phase (fyf 10° [40,44]. 

However, the NP effects could also enter through the B^ —>■ J/ipK^ amplitude, where 
EW penguin topologies, which have a sizeable impact on this decay [45], offer a partic- 
ularly interesting scenario. The S-factory data for B tttt, ttK modes may actually 
indicate a modified EW penguin sector with a large CP-violating NP phase through the 
results for mixing-induced CP violation in B^ t^^Kq, [16,46], thereby complementing 
the pattern of such CP asymmetries observed in other h ^ s penguin modes, where the 
B^ (pKs channel is an outstanding example [7]. The sign of the corresponding CP- 
violating NP phase would actually shift Aqp{B(1 J/ipKs) in the right direction [16,45]. 
The interesting feature of the -B(i(s) decays is that they are essentially unaf- 

fected by such a NP scenario as EW penguins contribute only in colour-suppressed form 
and play a minor role. Consequently, a difference between the values of (pd extracted 
from Bd J/ipKs and the B^^s) — ^ ^d{s)^d(s) system could reveal such effects. 

A similar comment applies to the determination of the B'^-B'^ mixing phase, where 
the "golden" strategy uses mixing-induced CP violation in the time-dependent angular 
distribution of the Bg J/ip[-^ K^K~] decay products [41,47]; penguin 

effects can be controlled with the help of Bd — > J/ipp^ [48]. This determination of 0^ 
could also be affected by CP-violating NP contributions entering through EW penguin 
topologies. On the other hand, the extraction discussed in Subsection 14.31 is essentially 
unaffected, so that a difference between the two results could again signal such a kind 
of physics beyond the SM. Moreover, also a simultaneous analysis of the f/-spin-related 
Bs{d) J/ipKo, decays should be performed [2]. In analogy to the discussion given 
above, the (small) penguin effects in the determination of 0^ from Bd J/ipK'g, can 
then be controlled, and 0^ could be extracted from the h ^ d channel B^ J/ipK^, 
again with a sensitivity to a modified CP-violating EW penguin sector. 

As was noted in Ref. [49], the analysis of the Bd(s) ^d(s)'^d{s) decays can also 
straightforwardly be applied to the Bd(s) K^K^ system. Following these lines, the 
penguin effects in the determination of sin0<i from the h ^ s penguin decay 5° — > K^K^ 
can be included through its B^ K^K^ partner [50]|f| this is also the case for the 
corresponding Bd{s) x*^K*^ decays [48,51]. Again in these transitions, EW penguin 
have a very small impact. Should the interesting pattern in the mixing-induced CP 
asymmetries of B^ — >■ tt^/Ts, B^ — > (pKg and similar modes originate from a modified 
EW penguin sector, we would again not see it in the Bd[s) ^ K^*^°K^*^'^ system. 

3Here -> K°K° and 5° ^ K°K° take the roles of B", D+D- and 5° ^ D+D-, respectively. 
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6 Conclusions 



The CP violation in D^D^ offers another interesting probe for the exploration of 

the Kobayashi-Maskawa mechanism of CP violation. In these studies, the penguin effects 
have to be controlled, which can be done with the help of the [/-spin-related DfD~ 
channel. Motivated by the recent data from the B factories and the Tevatron, as well 
as the quickly approaching start of the LHC, we have investigated the allowed region in 
the space of the mixing-induced and direct CP violation of the S° D^D^ decay, with 
useful results to monitor the future improvement of the experimental picture, and have 
performed theoretical estimates of the relevant hadronic parameters and observables. 

We then discussed the extraction of CP-violating phases, where we may either use 
(pd as an input to determine 7, or use 7 to extract (p^- Concerning the former option, 
the current data point towards an unstable situation for the extraction of 7, while the 
strong phase 6 could be well determined. It appears therefore more interesting to extract 
the B^-B^ mixing phase from the CP asymmetries of B^ D~^D~, also since precision 
measurements of 7 will be available from the LHCb experiment through other strategies. 
We have provided the formalism to include the penguin effects, and have illustrated its 
practical implementation. In the case of the CP asymmetries of the B^ D^Dj decay, 
the penguin effects are doubly Cabibbo-suppressed and play therefore a significantly less 
pronounced role. However, they can also be taken into account with the help of the 
B^ — > D^D~ decay, allowing then a precision measurement of the B^-B^ mixing phase 
from the mixing- induced CP violation in B^ —>■ D^D~. 

An interesting feature of these determinations is the fact that they are insensitive to 
CP- violating NP contributions entering through the EW penguin sector. In this respect, 
they are complementary to the well-known standard strategies. The determinations of 
the B^-B^g mixing phase through the Bs{d) Dt{d)-^^{d) system on the one hand and 
Bs J/i'(p, Bd J/ipp^ on the other hand are particularly promising, and the studies 
of LHCb in this direction should be further pursued to fully exploit the physics potential 
of these decays. 
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